Testing Partonic Charge Symmetry at a High-Energy Electron Collider 
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We examine the possibility that one could measure partonic charge symmetry violation (CSV) 
by comparing neutrino or antineutrino production through charged-current reactions induced by 
electrons or positrons at a possible electron collider at the LHC. We calculate the magnitude of 
CSV that might be expected at such a facility. We show that this is likely to be a several percent 
effect, substantially larger than the typical CSV effects expected for partonic reactions. 
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Charge symmetry is a very specific operation 
involving isospin, which leads to the interchange 
of protons and neutrons, or equivalently the inter- 
change of up and down quarks. The charge sym- 
metry operator Pes corresponds to a rotation of 
180° about the 2 axis in isospin space, such that 

Pes = , 
Pcs\u) = -\d); P cs \d) = |«). (1) 

It is of particular importance because at low ener- 
gies, where it has been studied extensively, charge 
symmetry is a far better symmetry than isospin in 
general, typically being respected to better than 
1% [U [2]. It is therefore natural to assume that 
charge symmetry is also valid at the partonic level 
and, indeed, almost all analyses of parton distri- 
bution functions (PDFs) assume charge symmetry, 
whether the assumption is stated or not. The im- 
portance of charge symmetry violation in PDFs 
within the context of tests of the Standard Model 
has recently been of considerable interest [21 0] . 

To date, no violation of charge symmetry has 
been observed at the partonic level, although the 
one global analysis that did allow for CSV did 
find a preferred solution with a non-zero effect 
- albeit with very large errors [5]- The current 
upper limits are consistent with the validity of 
partonic charge symmetry in the range 5-10% [BJ. 
Theoretical models generally produce estimates of 
charge symmetry violation (CSV) in PDFs which 
for many observables give effects at roughly the 
1% level [BJ [7]. This presents a significant chal- 
lenge for experimentalists, first to observe effects 
of this magnitude and then to isolate the signal 
from competing effects of similar size. 

A new facility has recently been proposed that 
would collide electrons or positrons from an elec- 



tron accelerator with protons or deuterons from 
the LHC (8] . In this paper we will show that such 
a facility (given the name LHeC) has the poten- 
tial to produce charge symmetry violating effects 
which are considerably larger than those expected 
with other facilities. We will review the effect in 
question, show the results of theoretical calcula- 
tions for the proposed CSV effects, and discuss 
why they ought to be expected to be relatively 
large at energies accessible to an electron-ion col- 
lider. 

The reactions of interest are the charged cur- 
rent (CC) cross sections for electron and positron 
deep inelastic scattering at energies in the range 
50-100 GeV on protons and deuterons at LHC 
energies, i.e., several TeV. These are important 
because they directly and unambiguously probe 
the flavor structure of the proton PDFs in the va- 
lence region. Consider the deep inelastic reaction 
(e~, u e ). A high-energy electron incident on a pro- 
ton produces a neutrino. The process results from 
a W~ which is absorbed on quarks from the pro- 
ton, as shown schematically in Fig. [TJ the final 
hadronic state is not observed. The signature for 
this process is disappearance of the electron, to- 
gether with very large deposition of energy in the 
hadronic sector. 

The i<2 structure function for the CC reaction 
on a proton has the form 

fY~ p (x) = 2x{u(x) + c(x) + d(x) + s(x)]. (2) 

These reactions will occur at extremely high ener- 
gies and very large Q 2 . Therefore, any corrections 
to the F 2 structure functions in Eqs. Q and ^ 
arising from quark mixing matrices, quark masses 
or higher-twist effects should be completely negli- 
gible. 
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FIG. 1: Schematic picture of charged-current neutrino 
production in DIS induced by an electron on a proton. 



We can also consider the corresponding reaction 
for positrons on protons, (e + ,V e ). This reaction 
involves the absorption of a W + on the proton, 
with the resulting F 2 structure function 

F 2 w+P (x) = 2x[u(x) + c(x) + d{x) + s(x)}. (3) 

We can straightforwardly calculate the F 2 struc- 
ture functions (per nucleon) on the deuteron, 

F 2 w ~ D (x) = x[u+(x)+d+(x) + 2c(x) 

+ 2s(x) - Sd(x) - 6u(x)] ; 

F 2 w+D (x) = x[u + (x) + d+(x) + 2c{x) 

+ 2s{x) - 6d(x) - 5u(x)} . (4) 

In Eq. Q we introduce combinations of quark 
parton distribution functions (PDFs) that are even 
or odd under charge conjugation, and the CSV 
PDFs 

q ± {x) = x[q(x)±q(x)] ; 
Su(x) = u p {x)-d n {x) ; 
Sd{x) = d p {x)-u n (x) . (5) 

There are analogous relations to Eq. ([5| for the 
antiquark CSV PDFs. For the remainder of this 
letter we assume that c(x) — c(x). The distribu- 
tions q~{x) 1 which involve the differences between 
quark and antiquark PDFs (alternatively, they are 
the C-odd combinations of quark distributions), 
are the valence parton distributions for a given 
quark flavor. 

Now define the following quantity, 



F, 



{x)+F 2 w+p {x) 



The quantity R~(x) is given by the difference in 
the F 2 structure functions per nucleon for electron- 
deuteron and positron-deuteron CC reactions, di- 
vided by the average F 2 structure function for CC 



reactions on protons initiated by electrons and by 
positrons. 

Using Eqs. §, <§ and @ we can straightfor- 
wardly show that the quantity Rr(x) in Q has 
the form 



R-(x) = 



x[— 2s (x) + Su (x) — Sd (x)] 
x[u+(x) + d+(x) + s+(x) + 2c(x)] 



(7) 

Thus R (x) is proportional to the valence quark 
CSV parton distributions plus the strange quark 
asymmetry (the difference between the strange 
and antistrange PDFs). Insofar as the strange 
quark asymmetry exists, it should be large only 
at quite small Bjorken x < 0.1, while theoretical 
estimates of the valence CSV parton distributions 
[! Eg suggest that for Q 2 - 10 GeV 2 they peak 
at values x ~ 0.4. 

Consider a hypothetical collider with 50 GeV 
electrons or positrons colliding with protons and 
deuterons of energy roughly 7 TeV. This would 
be similar to the possibilities if an electron col- 
lider were built at the LHC. Consider charged- 
current reactions at such a facility with Q 2 = 10 5 
GeV 2 . We know of two different mechanisms for 
charge symmetry violation in parton distribution 
functions. The first arises from the radiation of a 
photon by a quark. Such contributions are shown 
schematically in Fig. [2j they were first calculated 
by the MRST group [U and Gluck et al. [12]. 
These QED corrections are analogous to the cou- 
pling of gluons to quarks, except that photons do 
not have the self-coupling terms possessed by glu- 
ons. Inclusion of these 'QED splitting' terms will 
produce charge symmetry violation in parton dis- 
tribution functions because of the electromagnetic 
(EM) coupling due to the different charges of up 
and down quarks. 
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FIG. 2: Schematic picture of quarks coupling to pho- 
tons. This gives the origin of QED splitting that pro- 
duces CSV effects in parton distribution functions. 

The behavior of parton distributions with in- 
creasing Q 2 is given by the DGLAP evolution 
equations [TBT - TTB] . We expand the DGLAP evolu- 
tion equations to lowest order in both the strong 
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coupling a s and the electromagnetic coupling a, 
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d In /i 2 
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e 2 p lq ® ft 



(8) 



In Eq. ds]), qi(x,/i 2 ) is the parton distribution for 
a given flavor i, g(x,^L 2 ) is the gluon distribution 
and 7(2:, /^ 2 ) is a "photon parton distribution" [TT] . 
In these equations the convolution is defined as 



1 ^P(p)ff(-,M 2 ) . 

2/ y 



(9) 



and the splitting functions are given by 



P 99 (y) = 
P 9 7(y) = 



c F ' 
P w (y). 



P 79 (y) = 



Pqqiy) . 

c F ■ 



(10) 



We have analogous equations to Eq. pj) for an- 
tiquaries. Taking the valence combinations from 
Eq. jjj), we obtain for up and down valence quarks 



du~(x, n 2 ) 
d In ^i 2 

dd~(x, fj, 2 ) 
din jj? 



^P qq ® ti- 
^P M ®d- 



2a ~ 



(ii) 



For the valence CSV parton distributions, since 
Su~(x) = Up(x) — d~(x), from Eq. (fill) we obtain 
the evolution equations for the valenceCSV PDFs, 
to lowest order in q s and a, 



d[5u (x,fi 2 )] 

<91n/i 2 
d[5d-(x,fi 2 )} 

91n/i 2 
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(el ~ e 2 d )P q , 



- — (e 2 



el)P, 



m 



dr .(12) 



Eq. describes how the valence quarks evolve 
with Q^and Eq. (12 1 shows how the valence CSV 



distributions evolve, to lowest order in both a s and 
a. With increasing Q 2 , partons radiate gluons and 
photons which carry off momentum. Since the to- 
tal momentum fraction carried by quarks is given 
by the second moment of the parton distributions, 
as Q 2 increases the parton distribution functions 
will shift towards progressively smaller x values. 



Comparison of Eqs. (Ill and ( 12 1 shows that the 



radiation from valence quarks will be greater than 
that from the valence CSV distributions. This oc- 
curs because to lowest order in a s and a, valence 
quark evolution contains contributions from both 
gluon and photon radiation, whereas the valence 
CSV distribution has only a term from photon ra- 
diation. This suggests that with increasing Q 2 the 
valence parton distributions would experience a 
larger shift to low x than will the valence CSV 
distributions. We note that the quantity R~(x) 
defined in Eq. ^ is proportional to the ratio of 
valence CSV distributions to valence PDFs, at a 
given x value. If the CSV valence distributions 
are becoming larger relative to the valence PDFs 
at large Q 2 , then we expect the quantity R~(x) to 
grow as Q 2 increases; specifically we would expect 
the ratio to increase logarithmically with Q 2 . 

Eq. (12), the QCD evolution equations for 
the valence CSV parton distributions, have been 
solved by Gliick et at 12 and also by the MRST 
group [TT]; the two groups made slightly differ- 
ent approximations for the initial conditions. We 
stress that while the effect of photon radiation is 
clear, it is far less obvious that the boundary con- 
ditions imposed on the calculations are appropri- 
ate. That is, we know of no rigorous proof that a 
low scale, typical of quark models, is the appropri- 
ate place to set the effect to zero. In the absence of 
a compelling theoretical derivation it is extremely 
helpful to be able to test the idea experimentally. 

The MRST group [TT] did attempt an experi- 
mental test of this method. Including QED radi- 
ation in the DGLAP equations introduces a 'pho- 
ton parton distribution' j(x,Q 2 ), which appears 
in Eq. @. The MRST group attempted to iden- 
tify this quantity in the process ep — > ejX where 
the final state e and 7 are produced with equal 
and opposite large transverse momentum. This 
process has been measured by the ZEUS Collabo- 
ration in ep collisions at s/s — 300 and 318 GeV 
|16j . The observed cross sections were in reason- 
able agreement with the MRST calculations but 
disagreed with calculations done using the Monte 
Carlo simulations PYTHIA [JH and HERWIG 
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[IB]. It would be useful to have other experimen- 
tal tests of this method for including radiation of 
photons by partons, and the experiment suggested 
here could provide additional confirmation of this 
method. 

A second source of valence parton CSV arises 
naturally from the mass difference between the u 
and d quarks and may be calculated within light 
cone quark models. In such models the valence 
quark distribution can be expressed as [TOH2"T] 



g v (x,/x 2 ) = \(X\ 



1 + 7 U 7 



0„,3 



x S(M(l-x)-p+) . 



(13) 



Eq. ( 13 1 denotes the process where a valence quark 



is removed from a nucleon \N), and the result is 
summed over all final states \X). The quantity p\ 
is the energy of the state following removal of a va- 
lence quark with momentum k. The quantity [i 2 
represents the starting value for the Q 2 evolution 



of the parton distribution. Eq. (131 is formally ex- 



act and provides a natural starting point for calcu- 
lations which preserve the correct support of the 
PDFs. 

Model quark wavefunctions are found to be 
nearly invariant under the small mass changes typ- 
ical of CSV [TO], so we concentrate on the breaking 
of partonic charge symmetry associated with en- 
ergy shifts resulting from the u and d quark mass 
differences. In particular, we consider the effect of 
the n — p mass difference 6M = M n — M p = 1.3 
MeV, as well as the difference in diquark masses 
arising from the current quark mass difference be- 
tween up and down quarks. We define the quantity 



6m = m dd - m v 



(14) 



for which we have a robust estimate Sfh ~ 4 MeV 
|22j . We determine CSV valence PDFs by calcu- 
lating the variation of quark model parton distri- 



butions from Eq. (131 with respect to these quan- 
tities, i.e., 



dq v dq v 



(15) 



From Eq. ( 15 ) the valence charge symmetry vi- 



olating parton distributions are obtained by tak- 
ing variations with respect to diquark and nucleon 
masses on valence parton distributions from quark 
models. The resulting PDFs account for quark 
and nucleon mass differences that lead to CSV ef- 
fects. 



Valence CSV parton distributions arising from 
quark mass difference effects were calculated by 
Rodionov et al. [TO]. They used bag model wave- 
functions, including the effect of quark mass differ- 
ences on the quark wave functions as well as on the 



di-quark and nucleon masses, using Eq. (13). The 



Rodionov calculation preserved the correct sup- 
port and included the effect of transverse momen- 
tum in the proton. As in any quark model calcula- 
tion, the resulting leading twist PDFs are appro- 
priate to a relatively low momentum scale (where 
most of the momentum of the nucleon is carried 
by valence quarks [2T1EI]). ^ n or der to compare 
with experimental data these PDFs arc typically 
evolved up to Q 2 = 10 GeV 2 . We subsequently 
evolved these parton distributions to the higher 
Q 2 values appropriate to an electron-ion collider 
such as the LHeC. 

An alternate theoretical approach was due to 
Sather [9], who investigated the expression for va- 
lence parton CSV distributions in a static quark 
picture. In such models the correct support is 
no longer guaranteed. In addition, Sather ne- 
glected transverse quark momentum. By apply- 



ing Eq. (15) to Eq. (13) within this approxima- 



tion scheme, Sather obtained an analytic approx- 
imation relating valence quark CSV to deriva- 
tives of the valence PDFs. The analytic approx- 
imation of Sather is appropriate only at Q 2 val- 
ues appropriate for quark model calculations, i.e., 
Q 2 ~ 0.25-0.5 GeV 2 . 

We used the Sather prescription, differentiating 
valence parton distribution functions to obtain va- 
lence CSV PDFs. For this purpose we used the 
MRST2001 parton distributions [H] at the start- 
ing scale, Ql = 1 GeV 2 . This is slightly too large 
a value of Q 2 for the validity of Sather's analytic 
approximation, but the resulting errors should be 
small. We then inserted the resulting CSV PDFs 
into the DGLAP evolution equations and evolved 
to the Q 2 appropriate for the electron collider ex- 
periments. The results were similar to those ob- 
tained using the CSV distributions of Rodionov et 
al. . 

Since the CSV effects arising from QED split- 
ting effects and from quark mass differences are 
nearly independent, we have added the two effects 
to produce a "net" CSV effect. 

The quark model estimates of valence parton 
CSV can be compared with a recent lattice calcu- 
lation of valence charge symmetry violating parton 
distributions [25] • The lattice calculation provides 



4 



striking confirmation of the quark model results 
for parton CSV. The lattice calculation was car- 
ried out by considering small deviations from the 
SU(3) flavor-symmetric point where the strange 
and light quarks masses are all equal. In that way 
they could estimate the effects of quark mass dif- 
ferences on parton distributions. The lattice re- 
sults gave estimates for the second moment of the 
valence CSV distributions 

SU+ = [ x5u + (x)dx = -0.0023(6) , 
Jo 

5D + = / x5d + (x) dx = +0.0020(3) . (16) 
Jo 

The lattice results were obtained at a momentum 
scale Q 2 = 4 GeV 2 . Note that the results are 
appropriate for the C-even combination of quarks 
rather than the desired C-odd combination for va- 
lence quarks, because the lattice calculations are 
sensitive to the C-even combination. By com- 
parison, the quark model valence quark calcula- 
tions at a similar scale obtained 6U~ = —0.0014 
and SD~ = +0.0015 [TU1 US]- Note that the lat- 
tice result agrees with the quark model results in 
both the sign and relative magnitude of the second 
moment of the valence CSV distributions. The 
lattice results are 30-50% larger than the quark 
model values. The differences may result from 
the inclusion of singlet contributions in the lat- 
tice calculations. The lattice results are also in 
good agreement with the best value obtained for 
valence quark charge symmetry violation in a phe- 
nomenological global fit to high energy data by 
the MRST group [5]. However, the uncertainties 
in the lattice calculation are considerably smaller 
than those from the global fit. 

There is one final term that enters into the quan- 
tity R~(x) of Eq. Q, namely the strange quark 
momentum asymmetry [27, 28 

xs~(x) = x[s(x) — s(x)] . (17) 

Strange (antistrange) parton distributions can be 
measured through opposite-sign dimuon produc- 
tion initiated by neutrinos (antineutrinos). A neu- 
trino undergoes a charged-current reaction, pro- 
ducing a fjT and a W + , which is absorbed on an s 
quark producing a charm quark. The charm quark 
subsequently undergoes a semileptonic decay pro- 
ducing a fi + and an s quark. The cross section for 
this process is proportional to the strange quark 
distribution. The corresponding reaction initiated 
by an antineutrino measures the antistrange PDF. 



Dimuon cross sections have been measured by 
the CCFR [H] and NuTeV [30] experiments. 
From these reactions one can extract the quan- 
tity xs~(x). These analyses have been undertaken 
by five groups: CTEQ [31]; Mason et al. [3J; 
the NNPDF Collaboration [33] ; MSTW [33]; and 
Alekhin, Kulagin and Petti [35]- We used the re- 
sults of the analysis of Mason et al. of the NuTeV 
neutrino reactions [32] . We made an analytic fit to 
the best-fit result from Mason et al. corresponding 
to Q 2 = 16 GeV 2 . The fit had the form 

xs~ (x) = Ax b exp(-cx)(x - 0.004) . (18) 

The resulting strange quark asymmetry was in- 
serted into the DGLAP evolution equation and 
evolved to high Q 2 . 
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FIG. 3: [color online] Parton distributions that oc- 
cur in the numerator of Eq. fl7J. Solid curve: 
xSu~ (x); long-dashed curve: xSd~(x); short-dashed 
curve: xs~(x). The PDFs have been evolved to 
Q 2 = 10 5 GeV 2 . 

The parton distribution functions that occur in 
the numerator of Eq. Q are plotted in Fig. [3] 
The solid curve is xdu~(x), the long-dashed curve 
is x5d~(x) and the short-dashed curve is xs~(x). 
As one might expect, the valence CSV distribu- 
tions peak at a relatively large value x ~ 0.2 
while the strange quark asymmetry peaks at an 
extremely small x value. Note that due to valence 
quark normalization, all of these quantities must 
have zero first moment, i.e., (q(x)) — 0, where 
q = [6u~ , 5d~ , s~]. The strange quark asymme- 
try has zero first moment because the proton has 
no net strangeness; the valence CSV distributions 
must have zero first moment because otherwise 
this would change the total number of valence 
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quarks in the neutron. So each of these curves 
crosses zero at a small value of x (not shown in 
Fig.§. 

Another notable point is that the signs of these 
quantities are such that (for values of x above the 
crossover point for all of the parton distributions) 
all three contributions should add together in the 
numerator of Eq. Q . Fig. [4] shows the expected 
value of R~ (x) vs. Bjorken x. The solid curve in 
Fig. [4] includes only the QED splitting contribu- 
tion to partonic CSV. The long-dashed curve in- 
cludes both QED splitting and quark mass contri- 
butions to valence quark CSV. The short-dashed 
curve is the result including all three terms in the 
numerator of Eq. Q , including also the contribu- 
tion from strange quark asymmetry. 
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FIG. 4: [color online] Contributions to the quantity 
R~ (x) vs. x from Eq. ([7]), where the PDFs are evolved 
to Q 2 = 10 s GeV 2 . Solid curve: contribution from 
QED splitting parton CSV term only; long-dashed 
curve: includes contribution also from quark mass 
CSV term; short-dashed curve: contribution from all 
terms including strange quark asymmetry. 

We see that for large x > 0.2 the strange quark 
contribution is essentially negligible. The pre- 
dicted values of R~(x) are large; for x = 0.6 the 
calculated ratio is greater than 6%. This is quite 
a sizeable result for partonic CSV terms, which 
for most observables yield effects at the 1% level 
or smaller [5J. This confirms our argument that, 
while both the valence quark and valence CSV dis- 
tributions shift to lower x values with increasing 
Q 2 , the CSV distributions experience a smaller 
shift (because to lowest order the CSV valence dis- 
tributions only radiate photons while the valence 
parton PDFs radiate both gluons and photons), 
and thus for a given x the ratio of valence CSV dis- 



tributions to valence PDFs should increase slowly 
with Q 2 . Our best theoretical estimate of the ra- 
tio R~(x) from Eq. Q at large x values is pre- 
dicted to be rather large, of the order of several 
percent. For reasonably large values x > 0.1, the 
ratio R~(x) is composed of relatively equal con- 
tributions from valence parton CSV effects arising 
from quark mass differences and from QED radia- 
tion. Thus the quantitative values obtained for the 
ratio R~ (x) can provide a further check on the as- 
sumptions made in determining charge symmetry 
violation arising from QED radiation. 

In conclusion, a high energy electron/positron 
collider whose beams interact with deuteron 
beams from the LHC may produce the most 
promising observable with which to search for par- 
tonic charge symmetry violating effects. 
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